Quantum fluctuation and geometrical frustration effects on electric polarization 
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We examine theoretically a possibility of ferroelectricity caused by electronic charge order without inversion 
symmetry, motivated by layered iron oxides. Quantum electronic models in a paired-triangular lattice are ana- 
lyzed by utilizing the variational Monte Carlo simulation. Our calculation demonstrates that combined effects 
of electron transfer between the layers, corresponding to quantum fluctuation between the potential minima, 
and geometrical frustration promote appearance of an electric polarization. Present results are in contrast to the 
conventional manner of quantum fluctuation in ferroelectricity. 



PACS numbers: 77.80.-e, 71.10.-w, 75.80.+q 

In solid state textbook, it has been explained that ferroelec- 
tric transition in the displacive-type ferroelectricity is caused 
by a balance of the long-range Coulomb interaction and the 
short-range repulsion between ions. In addition to these clas- 
sical explanation, recently developed ab-initio band calcu- 
lations reveal that the electron covalency between cations 
and anions have considerable contributions in the ferroelec- 
tric transition. As for the order-disorder type ferroelectricity, 
the classical dipole-dipole interaction aligns the permanent 
dipoles and brings about the ferroelectric transition. Quantum 
effect plays irnportant roles on the hydrogen-bond type ferro- 
electricity |[ll 0] • The quantum proton tunneling between the 
two potential minima suppresses the transition, and arouses 
the quantum phase transition between the ferroelectric and 
quantum paraelectric phases. 

There is another class of ferroelectricity where an elec- 
tric dipole and a ferroelectric transition are attributed to the 
electronic origin. This is termed an electronic ferroelectricty. 
For this type of the ferroelectricity, several experiments in 
transition-metal oxides IHHEl] and organic salts 0,0,1111 have 
been reported, and some exotic theoretical ideas have been 
proposed JoKl^l- Most of the systems belong to the charge- 
ordered insulator with the so-called 1/4 filling, and a charge 
alignment without the inversion symmetry is responsible for 
the electric dipole moments. Here we give a simple exam- 
ple in Fig. [T] a charge order (CO) model in a pair of chains 
or layers. In the CO with the three-fold or odd-fold periodic- 
ity, a pair of the charge-neutral sites, represented by the grey 
circles, is sandwiched between the checker-board type charge 
configuration. Quantum charge fluctuations are caused by the 
electron transfer between the two potential minima. When the 
electrons are localized in one side due to some cooperative 
interactions, the macroscopic electric dipole appears. 

Instead of extensive material syntheses and experimen- 
tal examinations in the electronic ferroelectricity, the micro- 
scopic mechanisms have not been clarified yet. The long- 
range Coulomb interaction may promote an electronic charge 
ordering associated with electric polarization. In comparison 
with other types of feiToelectricity, quantum fluctuation due 
to electron motion should be crucially important, because the 
electron mass is much smaller than the proton and ion masses. 
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FIG. 1: (Color online) A three-fold charge order structure and the 
electric polarization. Black, white and gray circles represent charge 
rich, poor and neutral sites, respectively. 

Large fluctuation causes the diffusive nature in the dielectric 
constant and the opticalspectra observed commonly in this 
class of material Jj, 0, [Till . It is expected naively that such 
fluctuation prevents a system from the long-range dielectric 
order, on the analogy of the quantum phase transitions in the 
hydrogen-bond type ferroelectricity and SrTiOa. 

In this Letter, we demonstrate, against such conventional 
expectation, that, in a frustrated geometry, a quantum charge 
fluctuation promotes the electric polarization. The present 
study is motivated by an electronic ferroelectric material, 
the layered iron oxide i?Fe204 (i?=Lu, Yb, Y). The crystal 
structure is made up of a stacking of the paired-triangle lat- 
tices 11211 . as shown in Fig.|4ta). The nominal valence of a Fe 
ion is 2.5h- implying a coexisting of an equal amount of Fe^+ 
and Fe'^^. The electric polarization measured by the pyroelec- 
tric current is enhanced around the charge ordering tempera- 
ture of the three-fold 10,[l3l[l3l- A charge disproportionation 
in the paired-triangle layers is believed to be responsible for 
the electric dipole moment, as expected in Fig. [T] Here we 
concentrate on a paired-triangle lattice, and take up a possi- 
bility of the electronic ferroelectricity. In particular, as well 
as the quantum charge fluctuation, we focus on roles of the 
geometrical frustration which is often seen in the electronic 
ferroelectric materials iH, 01- To explore this issue, the 
three electronic models are analyzed by the numerical method 
of the variational Monte-Carlo (VMC) simulation. We show 
that the electron transfer between the layers, corresponding 
to the quantum fluctuation between the potential minima, sta- 
bilizes the electric polarization, when the system is strongly 
frustrated. 

We start from a possibility of the electronic ferroelectric- 
ity in a simple spin-less fermion model where the long-range 
Coulomb interactions and the electron transfers are taken into 
account. We consider a pair of the triangular layers where 
two triangular lattices are regularly stacked along the c axis as 
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shown in Fig.|2a); a site in the upper layer is just upon a site 
in the lower one. Effects of the realistic crystal structure of 
i?Fe204 and spin degree of freedom will be considered later. 
The Hamiltonian is given as Tivt ='Ht + Ti-v with 



nt = - t 



<ij>m 



H.c. 



(1) 



and 



Hv = nmirimj + V'y^ 

+ VlV" nuiUii, 



{ij)m 



(2) 



where c\^^ is the creation operator for a spin-less fermion at 
site i on the upper layer [m 
and 7i„ii = cl^jC 
a two-dimensional coordinate in a triangle lattice. To exam- 
ine roles of frustration, we consider an anisotropic triangular 
lattice where bonds along the [110] direction are inequivalent 
to those along [100] and [010] [see Fig.l^a)]. Introduced V, 
V' and V± are the Coulomb interactions between the nearest- 
neighboring (NN) sites along the [100] and [010] directions, 
the [110] one, and the [001] one, respectively. Symbols J2{ij) 
and imply summations for the NN pairs. The three 

transfer integrals, t, t' and t±, are defined in the same way. 
Ratios V' /V and t' /t express measure of frustration. 

Our optimization VMC method can perform unbiased cal- 
culation where both the frustration and quantum effects are 
treated accurately lIlSll . We adopt 2-5 x 10^ samples in most of 
the simulations where an error of the energy expectation value 
is of the order of 10~^t. Systems used are of Ng = L x L x 
2(= 2iV) (L < 12) sites with the periodic- and antiperiodic- 
boundary conditions, and the particle number is = Ns/2. 
A variational wave function is defined as ^I* = 7^$, with the 
one-body Hartree-Fock (HF) part $ and the many-body cor- 
relation factor v. We define V = Ylij (1 ~ Vij^iUj) where 
Vij takes one of the three variational parameters, v, v' and v±, 
when the NN i j pairs are along the [100] and [010] direc- 
tions, the [110] one and the [001] one, respectively. As for $, 
we assume the HF wave functions for the following four-types 
of the CO structures; (i) a two-fold CO along the [110] direc- 
tion termed CO1/2-I, (ii) a two-fold CO along [100] termed 
CO1/2-II, (iii) a four-fold CO along [110] termed CO1/4, and 
(iv) a three-fold CO along [110] termed CO1/3 [see Fig.^a) 
and the inset of Fig. |2b)]. In each CO, a magnitude of the 
order parameter is introduced as a variational parameter A in 
3". In the CO1/2-I, COi/2-n and CO1/4 structures, stackings 
of the CO patterns in the two triangular-layers are set to be 
out-of -phase. In CO1/3, the assumed pattern along the [110] 
direction is schematically ■•••o*»o*-- - (•••o»*o»*---)in 
the upper (lower) plane. Symbols •, o and ★ imply the charge 
densities of {n.mi — 1/2) = A, —A and zero, respectively. 

In Fig. 12 a), we present the CO phase diagram. We set up 
t = t' = t± and V ~ V±, and focus on roles of the frustra- 
tion, V' /V, and the electron transfer In the classical limit. 



u) or the lower one (m — I), 
is the number operator A subscript i takes 
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FIG. 2: (Color online) (a) A paired-triangular lattice adopted in the 
Hamiltonian TLvt- (b) Phase diagram for Tivt- Relations t — t' = 
t± and V = V± are assumed, (c) The energies for the several CO 
patterns at V'/V — 1. (d) The polarization correlation P. Cir- 
cles, squares and inverse triangles are for the results in CO1/3 at 
V'/V = 1, COi/2 at V'/V = 0.95 and CO1/4 at V'/V = 1.05, 
respectively. Open and filled symbols are for the results in L=6 and 
12, respectively. 

i.e. t = 0, the CO1/2-I (CO1/4) structure are realized in the 
regions ofV'/V < 1 (> 1). The COi/2-n structure is degen- 
erate with CO1/4. The CO1/3 structure of the present interest 
appears only at V'/V = 1 where all CO's are degenerate. 
By introducing t, CO1/3 appears in a finite parameter region. 
That is, COi/3 is stabilized by a combined effect of the quan- 
tum fluctuation and the geometrical frustration. 

The energy expectation Stot = N^^{Hvt) is presented as 
a function of t/V at V'/V = 1 in Fig. |2|c). It is clearly 
shown that a stabilization of CO1/3 is attributed to the i-linear 
dependence of i^tot, in contrast to the quadratic dependence 
in other CO's. A schematic pattern of CO1/3 is shown in 
Fig. [21 a). The energy gain is caused by the inter-layer elec- 
tron transfer, t± as follows. Let us focus on a pair of the gray 
circles in the upper and lower layers in Fig. [2a). This system 
corresponds to the symmetric double-wells aligned on a tri- 
angle lattice. Since these sites are surrounded by the in-plane 
three NN charge-rich and three poor sites, the Coulomb inter- 
actions are canceled out. Thus, the classical energy does not 
depend on the charge configurations at the sites represented 
by the gray circles in Fig. [2a). This degeneracy of the order 
of 2^/"^ is lifted independently at each site by the first order 
of t±. Even when the longer-range Coulomb interactions than 
V, V and V± are considered, this scenario survives in the 
case where t± is larger than those interactions. 

The electric-polarization correlation defined by P = 

1 /2 

with p = {3/N)Y,i{nui - nu) is presented in 
Fig. [2d). The possible maximum value of P in CO1/3 is 
one. The obtained value in CO1/3 is much higher than those 
in other CO's. However, with increasing the system size L, 
a magnitude of P in CO1/3 decreases considerably. We sup- 
pose that amplitude of P in the thermodynamic limit is much 
smaller than its maximum value. 

The above results imply that the quantum charge fluctu- 
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ation and the geometrical frustration promote the three-fold 
CO upon which the electric polarization has a chance to ap- 
pear. However, in the large system size, the amplitude of P 
in the polar state is rather small. When we focus on one of 
the paired-grey sites in Fig.|2|a), corresponding to a symmet- 
ric double-well potential, the bonding state is stabilized by 
0{t±), and a symmetric charge distribution is realized. The 
present numerical results imply that, within this model, the ef- 
fective interactions between the pairs are not enough to induce 
an asymmetric charge distribution. We expect that additional 
small factors, which are not taken into account so far, may 
stabilize the polar CO structure. We propose two plausible 
candidates, the spin degree of freedom and the realistic crys- 
tal structure in i?Fe204. 

We are modeling the spin degree of freedom on the basis of 
LuFe204 where two kind spins, S = 2 (dP) in Fe^+ and S = 
5/2 (cf) in Fe'^"'', coexist. An excess electron in a Fe^+ is itin- 
erant and is responsible for CO, and others are regarded as a 
locaUzed spin with S = 2 because of their small transfer inte- 
gral M 

12111 . The model is made up of the itinerant electrons 
coupled ferromagnetically with the localized spins through 
the on-site Hund coupling. This is a kind of the double- 
exchange model. We assume, for simplicity, the infinite-limit 
of the Hund coupling and the Ising-type localized spins de- 
noted by <j„ii (m = u, I) which takes ±1. The Hamiltonian is 
Tis = 'H[+ Hv + "Wafm where Hv is given in Eq. Q. The 
first term is the same form with Eq. ([U, but the all transfer 
integrals are replaced by the spin dependent forms. For ex- 



ample, 
where 



<ij>m ^mi'^rnj " ' S<ij>m '•'^Tni'^mj '-mi^mj 

_ = t for cTjni = cTmj and zero for ami ^mj- 

The last term T^afm is the antiferromagnetic interaction be- 
tween the localized spins, Hafm = Js{Y.m<ij> '^mi'^mj + 
^{■ij)rn '^nii<yrnj+J2i '^uiCTii} with the positivc exchange con- 
stant Js- We analyze this model by using the VMC method. 
We consider all possible magnetic orders in the localized- 
spin sector characterized by the momenta (27r/3, 27r/3, 0) 
which corresponds to the observed neutron diffraction peak 
at (1/3 1/3 m) in LuFe204 S [13, ES [Hi ■ We restrict our 
calculations to the most frustrated case, i.e. V' /V = 1 and 
t'/t = 1. 

In Fig. [3 a), the energy and the polarization correlation in 
the optimized spin state are plotted as a function of N^/Ne- 
Here, A^kd is the number of the conduction electrons with 
spin up (down) and A^e = + -^i ■ The energy minima are 
seen at N^/N,, = 0.33 and 0.67 (NjNe = 0.33), where a 
large value of P {^^ 0.37) is realized. A magnitude of P is 
almost independence of the system size. That is, the polar 
state becomes robust by the spin ordering. 

This stabilization of the polar state is attributed to the so- 
called gate effect due to spin degree of freedom as follows. 
The optimized spin structure is presented in Fig. Ob). It is 
easily seen that, without the spin degree of freedom, the classi- 
cal energy is not changed by an electron motion from a site in 
the lower layer to the one just upon it. The energy gain in this 
process is proportional to t±, being similar to the previous sit- 
uation in the CO1/3. In contrast to the previous case, the spin 
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FIG. 3: (Color online) (a) The energy and the electric polarization 
correlation in the model with the spin degree of freedom, Tis- Re- 
lations t = t' = t± and V = V' = V± are assumed. Parameter 
values are chosen to be t/V = 0.2 and Js/t = 0.05. Open and 
filled symbols are for the results in L = 6 and 12, respectively, (b) 
A schematic view of the optimized spin structure. Short and long 
arrows represent spin directions of the conduction electrons and the 
localized spins, respectively. 

order suppresses the electron hopping, and keeps the electrons 
in the lower layer. This effect is seen, for example, for the 
down-spin electrons in the lower layer; since the number of 
the down localized-spins in the lower layer is twice of that 
in the upper one, the down-spin electrons are confined in the 
lower layer to gain the in-plane kinetic energy. In other word, 
the itinerant electrons are magnetically polarized by the local- 
ized spins. As a result, charge imbalance is induced between 
the layers, and the electric polarization appears. It is worth 
noting that the spontaneous spin- and charge-polarizations in 
the conduction electrons are not caused by an asymmetric po- 
tential due to the localized spins; the optimized localized-spin 
configuration is invariant under the space and time inversions 
plus the in-plane translation along [100]. 

Another candidate to realize the polar CO state is the realis- 
tic crystal structure in i?Fe204. Unlike in the models consid- 
ered so far, in i?Fe204, a position of a Fe ion in the upper layer 
is not just upon a Fe ion in the lower layer, but is upon an O ion 
surrounded by three in-plane Fe ions [see Fig.Ua)]. We mod- 
ify the previous Hamiltonian Hvt by considering this crys- 
tal structure. We introduce the long-range Coulomb interac- 
tions for the in-plane NN sites (KibNN), that for the inter-plane 
NN ones (Knn), and for the inter-plane NNN ones (VcNNn)- 
When the 1 /r-type Coulomb interaction is assumed, we have 
VcNN/VabNN = 1-2 and VcNNN/VabNN = 0.77 UM- As for 
the electron transfer term Ht, the transfer integral for the in- 
plane NN sites (tab), and that for the inter-plane NN ones (tc) 
are introduced. This model is analyzed by the VMC method. 
Instead of the CO1/3 introduced previously, we adopt the po- 
lar three-fold CO (polar CO1/3) and the six-fold one (COi/g) 
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FIG. 4: (Color online) (a) A paired-triangular lattices in i?Fe204 
and the polar CO1/3 structure. The lower panel shows a side view 
along the [120] direction, (b) Phase diagram for the model where the 
realistic crystal structure is taken into account. Relations tc/tab ~ 1 
and 'KNN/VabNN ~ 1 are assumed, (c) The electric polarization 
correlation P at 14NNN/VabNN = 0.5. Open and filled symbols are 
for the results in L = 6 and 12, respectively. 



which are degenerate in the HF level. The CO patterns along 
[100] is schematically ••••o»»o»-- - (•••o»oo»o-- - )in the 

upper (lower) plane for the polar CO1/3, and is ■ • ••o»»oo • • • 
(•■•o»oo»»-- - ) for COi/g, which is regarded as an antiferro- 
electric state. Symbols o and • imply the charge densities of 
{nmi — 1/2) = A and — A, respectively, with the variational 
parameter A. 

The phase diagram is shown in Fig. |4|b) in a plane of 
VcNNN/VabNN andtab/KbNN- Relations VcNN = KbNN and 
^ab — tc are assumed, although the qualitative features intro- 
duced below are robust around these relations. In the classical 
limit of tab = 0, VcNNN/14bNN = 0.5 is the fully frustrated 
point, where the polar CO1/3 is degenerate with other CO's. 
This is similar to the results in Fig. |2tb). By taking the elec- 
tron transfer into account, the polar CO1/3 appears around this 
point. Beyond tab/VabNN ~ 0.15, the non-polar COi/g over- 
comes the polar CO1/3. The electric polarization presented 
in Fig. Iljc) shows almost the maximum value in the polar 
CO1/3, and its value is independence of the system size L. 
On the other hand, a value of P in COi/g tends to vanish in 
the thermodynamic limit. 

Key interactions to stabilize the polar CO1/3 are the inter- 
layer charge transfer, t^, and the long-range Coulomb interac- 
tion, K:NNN, which is not disregarded in this crystal structure. 
As shown in Fig. Eta), T4nnn induces the interaction along 
the [120] direction between the bonds which are represented 
by the pairs of the black-white circles connected by the grey 
bars. These bonds correspond to the pairs of the grey circles 
in Fig. 12 a). As a result, asymmetric charge distributions in 
the bonds overcome the symmetric ones in the bonding states 
which are realized in the previous model of Hvt- Although 
the dipole moments are induced in the local bonds, there is 
still competition between the polar CO1/3 and the antiferro- 
electric COi/g which are in the different configurations of the 



dipole moments along [100]. This degeneracy is lifted by t±. 
Consider the charge fluctuation at the sites represented by the 
black (white) circles in the upper (lower) layer in Fig.Ua). In 
the inter-layer charge fluctuation by t±, the energy reduction 
is of the order of i^/FcNNN where the intermediate-state en- 
ergy VcNNN is the smallest Coulomb interaction in the present 
model of Hvt- Therefore, this kind fluctuation easily occurs 
and stabilizes the polar CO1/3. In other CO's, this fluctuation 
process is suppressed or prohibited. 

Through the present analyses, of the three-type models in 
a paired-triangular lattice we demonstrate that the quantum 
charge fluctuation due to the electron transfer promotes the 
electric polarization. In particular, the polar phase is stabilized 
remarkably due to the inter-layer charge fluctuation in the 
strongly frustrated condition. The spin degree of freedom of 
electron and the realistic crystal structure in i?Fe204 are also 
the key ingredients for the electric polarization. In compar- 
ison with the hydrogen-bond type ferroelectricity where the 
proton tunneling suppresses the long-range order, the charge 
fluctuation in the present electronic ferroelectricity plays a 
distinguish role. From the present scenario of the ferro- 
electricity driven by the quantum fluctuation, the experimen- 
tally observed diffusive dielectric constant jj] and the valence 
fluctuation sug gested by the optical and Mossbauer measure- 
ments il III22II are naturally understood. Here, a classical pic- 
ture, i.e. the polar charge order stabilized by the Coulomb in- 
teraction, is not applicable, and a fragile electric polarization 
is of prime essence. The present results open a new aspect of 
the quantum fluctuation in ferroelectricity. 

Authors would Uke to thank S. Mori, N. Ikeda 
H. Takashima and J. Nasu for their valuable discussions. This 
work was supported by JSPS KAKENHI, TOKUTEI from 
MEXT, and Grand challenges in next-generation integrated 
nanoscience. 



[1] R. Blinc et ai, J. Phys. Chem. Solids 13, 204 (1960). 

[2] P G. de Gennes et ciL, Sol. St. Comm. 1, 132 (1963). 

[3] K. Kato, and S. lida J. Phys. Soc. Jpn. 51, 1335 (1982). 

[4] N. Ikeda et at.. Nature 436, 1 136 (2005), and N. Ikeda et at., J. 
Phys. Soc. Jpn. 69, 1526 (2000). 

[5] Y. Tokunaga et al.. Nature Mat. 5, 937 (2006). 

[6] P Monceau, et al, Phys. Rev. Lett. 86, 4080 (2001). 

[7] K. Yamamoto et al, J. Phys. Soc. Jpn. 77, 074709 (2008). 

[8] S. Horiuchi and Y. Tokura, Nature Mat. 7, 375 (2008). 

[9] D. V. Efremov et al. Nature Mat. 3, 853 (2004)^ 

[10] J. van den Brink, and D. I. Kohmskii, arXiv:0803.2964l 
[11] X. S. Xu etal, Phys. Rev. Lett. 101, 227602 (2008). 
[12] N. Kimizuka et al, in Handbook on the Physics and Chemistry 
of Rare Earth edited by K. A. Gshneider Jr. and L. Eyring, 13, 
283 (Elsevier, 1990). 
[13] Y. Yamada et al, J. P hys. Soc. Jpn. 66, 3733 (1997). 
[14] A. M. Mulders et a/., [^av:0812.3094| 
[15] C. J. Umrigaref a/., Phys. Rev. Lett. 60, 1719 (1988). 
[16] J. Akimitsu et al, Sol. St. Comm. 32, 1065 (1979). 
[17] K. Shiratori et al. Proceedings of the 6th international confer- 
ence of Ferrites, 703 (1992). 



5 



[18] M. Angst etal, Phys. Rev. Lett. 101, 227601 (2008). 

[19] W. Wu, etai, Phys. Rev. Lett. 101, 137203 (2008). 

[20] A. Nagano et at., Phys. Rev. Lett. 99, 217202 (2007). 

[21] M. Naka etai, Phys. Rev. B 77, 224441 (2008). 



[22] S. Nakamura et al, J. Alloys Com. 275-277, 574 (1998). 

* Present address: Chiba Institute of Technology, Tsudanuma, 
Chiba 275-0016, Japan. 



